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We study transport in three dimensional Weyl semimetals with A'' isotropic Weyl nodes in the 
presence of Coulomb interactions or disorder at temperature T. In the interacting clean limit, 
we determine the conductivity by solving a quantum Boltzmann equation within a 'leading log' 
approximation and find it to be proportional to T, upto logarithmic factors arising from the flow of 
couplings. In the noninteracting disordered case, we compute the finite-frequency Kubo conductivity 
and show that it exhibits distinct behaviors for uj <^T and u) ^ T: in the former regime we recover 
the results of a previous analysis, of a finite conductivity and a Drude width that vanishes as NT^; 
in the latter, we find a conductivity that vanishes linearly with uj whose leading contribution as 
T — >■ is the same as that of the clean, non-interacting system a{uj,T = 0) = ^ . We compare 

our results to experimental data on Y2lr207 and also comment on the possible relevance to recent 
transport data on Eu2lr2 07. 



There has been a surge of recent activity studying 
Dirac excitations in two dimensional media, most fa- 
mously graphensi. A natural question is whether there 
are analogs in three dimensions, with a vanishing density 
of states at the chemical potential and linearly dispersing 
excitations. It has long been known that touchings be- 
tween a pair of non-degenerate bands are stable in three 
dimensions, and typically have linear dispersion. Near 
these, electronic excitations are described by an analog 
of the Weyl equation of particle physics, which describes 
two-component chiral fermions^r— . Hence these states 
have been dubbed Weyl semi-metals (WSMs)^. 

To remove a band touching (or Weyl node) one neces- 
sarily must connect to another node. This is in contrast 
with two dimensions: graphene's nodes can be gapped 
by different intranode perturbations that break inver- 
sion (I) or time reversal (7") symmetry. The enhanced 
protection in three dimensions is due to a topological 
property of the nodes - they are sources (monopoles) of 
Chern flux in the Brillouin zone (BZ). This momentum 
space topology is associated with several physical phe- 
nomena. In particular, it was recently realized^ that 
unusual surface states will result as a consequence of 
the band topology. These take the form of Fermi arcs 
that connect the projections of the nodes onto the sur- 
face BZ. Such topological properties are sharply defined 
as long as one can distinguish band touching associ- 
ated with opposite Chern flux. The presence of trans- 
lation symmetry, and hence conserved crystal momenta, 
is sufficient to protect these defining properties since the 
nodes are separated in the BZ. In principle one needs 
perfect crystalline order to define these phases; in prac- 
tice, smooth disorder that only weakly mixes nodes is 
expected to have little effect. Other manifestations of 
the band topology include an anomalous Hall effect^i 
that is tied to the momentum space displacement be- 
tween nodes, and magneto-resistance arising from Adler- 
Bell-Jackiw anomaly of Weyl fermions^i^. 

Physical realizations of WSMs require non-degenerate 
bands to touch; therefore spin degeneracy must be lifted 



(by either spin-orbit interactions or magnetic order) , and 
either 7" or I must be broken: otherwise, all bands 
would be doubly degenerate. We further require that 
the Fermi 'surface' consists exactly of the Weyl nodes. 
In 7"-breaking realizations where I is unbroken, a simple 
'parity criterion' applied to eight 7~-invariant momenta 
in the BZ can be used to diagnose the existence of Weyl 
nodes^. InS , certain pyrochlore iridates A2lr207 (A=Y or 
Eu), were proposed to be magnetically ordered WSMs, 
with N = 24 Weyl points, all at the Fermi energy^ 
reached similar conclusions but with N = 8. Alternate 
proposals include HgCr2Se4 in the ferromagnetic stateii 
and topological insulator-ferromagnet heterostructuresi, 
with N — 2, the minimum allowed. 

Motivated by the availability of transport data on the 
iridatesi^iiS,, we study the electrical conductivity of an 
idealized model of a WSM, with an even number N of 
isotropic Weyl nodes characterized by the same disper- 
sion, with N/2 nodes of each chirality as required by 
topology^iii. In the absence of impurities and inter- 
actions we expect the free fermion result, aQ^\uj) = 
N-nrr — - We demonstrate how this is modified in two 
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cases: 

(i) in clean undoped systems with Coulomb interactions, 
current is carried equally by counterpropagating elec- 
trons and holes and can be relaxed via interactions alone. 
Solving a quantum Boltzmann equation (QBE) we find 
a finite conductivity proportional to the temperature T 
(upto logarithmic factors), as expected of a quantum crit- 
ical system^^, where T is the sole energy scale. 
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and fine structure constant renormalized to the scale of 
the temperature T, and vp and ao — e'^/ehvp are the 
corresponding 'bare' values at the microscopic scale (See 

App. ig. 
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(ii) In the presence of impurities, power counting shows 
that white-noise disorder is an irrelevant perturbation, 
and a naive expectation is that the clean result ctq^'* is 
reproduced. However, the result is more interesting: by 
evaluating a standard Kubo formula, we find that the the 
finite-frequency conductivity exhibits different behaviors 
for uj <^T and uj ^ T: in the former regime we find in 
agreement witli^ a finite Drude-like response with a peak 
width vanishing as NT'^; in the latter, we recover a\^^ as 
the leading behavior, which is universal and independent 
of disorder. We also determine the manner in which the 
conductivity interpolates between these limits. 

Previous studies of 3D Dirac points have assumed 
Lorentz invariance^^ or worked at a topological phase 
transition between insulators^i. Although our work 
differs from both of these situations - instantaneous 
Coulomb interactions break Lorentz invariance, and we 
study a stable phase - there are sufficient parallels that 
a similar 'leading log' approximation suffices to solve the 
QBE. Coulomb interactions also lead to a finite dc con- 
ductivity in clean graphene - the 2D analog of a WSM 
- but the leading log approximation fails here and more 
analysis is neede d^^i^^ . 

The Model.- In a WSM, the electronic dispersion about 
a Weyl node is generically of the form i?wcyi = u ■ 
kl + X]a=i '^'^ ■ fecr^i where cr'^ are the Pauli matrices. 
The velocities satisfy v-^ ■ {v'^ x v^) ^ 0, and the Chern 
number ±f ('chirality') associated with the Weyl node 
is Sign (v^ • {v^ X v^)). For simplicity, we shall drop 
the term proportional to identity and assume isotropic 
dispersion; relaxing this assumption should only pro- 
duce small corrections. The Hamiltonian for a system 
of N identically dispersing Weyl nodes ('flavors') with 
Coulomb interactions and disorder may then be written 
a.s H = Ha + Hj + Ho, with (repeated indices summed) 

Ho=^Ha=^ / ^l^^ {XaVpk ■ a) 'lpk,a 
^ J kik2q 

Hd^ f J2i^l{r)U{r)Mr) (2) 

where 'ipk,a is a two-component spinor in the (pseudo)spin 
indices cr, a', a,b = 1 . . . index the flavors, vp is the 
Fermi velocity, which we set to unity, Xa = ±1 is the chi- 
rahty of the a*'* Weyl node, V{q) = ^p- describes the 
Coulomb interaction in a material with dielectric con- 
stant £, U{r) is a random potential with white-noise cor- 
relations {{U{r)U{r'))) = niVQS{r — r') where vq char- 
acterizes the strength of the individual impurities and rii 
their concentration, fi^ = J and we have written 

Hq assuming that the Fermi level is at the Weyl nodes, 
which is the only case studied in this Letter. Here and 
below we set h — kg — \e\ — f , and define /? = 1/T. 
Conductivity with Interactions.- Critical systems - such 
as graphene and the WSM at neutrality- are exceptions 



to the rule that disorder is essential for a finite conductiv- 
ity, since they support current-carrying states in which 
particles and holes transport charge with no net momen- 
tum by moving exactly opposite to each other. In con- 
trast to conventional finite-momentum charge transport, 
such deviations from equilibrium can relax in the pres- 
ence of interactions alone, leading to a finite conductivity. 

We study transport in an interacting WSM by solving 
a QBE for the thermal distribution function of quasipar- 
ticle states. In doing so, it is convenient to first calculate 
the current from a single node (but interacting with all 
the nodes) , before making the leap to the current carried 
by all N nodes. We focus on a node with fiavor a, which 
we take to have Xa = 1- The single-quasiparticle states 
are obtained by diagonalizingi/o: V/e,a C^Vfc.a = Jk,a, 
Ha UHaU^^ = Jk^'^pkjl^^'jkXa, and are labeled by 
their helicity A (the eigenvalue of a ■ p.) From now on 
we will suppress the index a. In general, operators cor- 
responding to various transport properties are not diag- 
onal in the helicity; diagonal contributions correspond 
to motion of particles and holes in the applied field and 
may be characterized by appropriate distribution func- 
tions f\{k,t) = (7fc;^7fcA), while the off-diagonal terms 
{l\fSi~\k) describe the motion of particle-hole pairs, 
which are finite-energy states. At particle-hole symmetry 
and for uj ^ T, contributions of the latter to transport 
are expected to be small, and we drop them forthwith. 
In this approximation, it is therefore sufficient to solve 
the QBE for quasiparticle distribution functions f\{k,t), 
subject to an external force F, 



(^^^+F-\/k^ fx{k,t)^~w[Mk,t)] (3) 



where w is the rate at which quasiparticles scatter out 
of the state (A, fc) at time t, and captures the effect of 
interactions. Our goal will be to determine the steady- 
state form of the non-equilibrium quasiparticle distribu- 
tion function. We will restrict ourselves to linear response 
in i^, i.e. we assume that the deviation of f\ from equi- 
librium is small. The result is a linear functional equation 
which may be recast as a variational problem. We solve 
the latter approximately by identifying 'leading log' con- 
tributions, which dominate the relaxation of the observ- 
able under consideration. As mentioned, we assume that 
the constants that enter the solution of ^ are renormal- 
ized to the energy scale of interest, namely T. 

Neglecting particle-hole pair contributions, the cur- 
rent is J{t) = -Jk{i^i(T^k)f -EA=±/fe-^^/A(fc'^)- 
For a weak applied electric field E{t), the deviation of 
f\{k,uj) = J dtf{k,t)e^'^* from the equilibrium distribu- 
tion function /°(fc) — (l -\- e'^^'^) ^, and hence the con- 
ductivity cr(cLi,T), can be parametrized^^ in terms of a 
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dimensionless, isotropic function g{k,uj): 

h{k,uj)^2^5{u:)fl{k) 

+\p^k.E{uj)[fl{k)f\{k)]g{k,^) (4) 

k 
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It therefore remains only to determine the function 
g(k,uj), to which we now turn. Inserting Q into 
and working to Hnear order in E, we find 



it(3L,gik,u) + l)flik)f_{k)k ^ C g{k,u)k 



where C 



is the colhsion operator, a Unear functional of g{k,uj)k 
given in [X] This is equivalent to the variational problem 
of extremizing the quadratic functionaU^""— 
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in which we have rescaled all momenta and frequencies 
by T. A key simplification, known as the 'leading log' 
approximation (LLA) stems from the power-law nature 
of the Coulomb interaction: as a result of this, logarith- 
mically divergent small-momentum scattering dominates 
C. We may write C = Co + Ci, which when thought of as 
linear functional of gk have eigenvalues of 0{a^ log a) 
and 0(0^), respectively. In the LLA we approximately 
optimize Q by choosing gk in the space spanned by eigen- 
states of Co; as shown in App. |^the choice g — k£,{uj) 
yields 
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Finally, we observe that the fiipped chirality of half the 
nodes is unimportant as they all give the same contri- 
bution to (t; thus using the result for S,{uj) in ([5]) and 
multiplying by N we find the result for N nodes^S , 



-i-fi^6.6 + Na'^ In a- 
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Note that in the case of graphene, the LLA fails because 
the log divergence stems from a phase space effect due 
to enhanced scattering of collinear particles, which can- 
not relax a current. Thus, the eigenstates of Co do not 
contribute to the relaxation, which therefore occurs only 
via subleading, noncoUinear scattering, i.e. Ci— li^. In 
3D, Co includes noncollinear and thus current-relaxing 
processes, so that the LLA analysis is sufficien t^^i^^ . 

In the dc limit, @ reduces to ([T]), which we may ratio- 
nalize in terms of the relaxation time approximation to 
the QBE. Here, we take 'w[fx{k,t)]] 



and use 



a combination of dimensional analysis and a straight- 
forward application of Fermi's golden rule to estimate 
the scattering rate as ~ Na^T. From this, we 
find the conductivity to be similar to ([1]), modulo log- 
arithms. This provides an estimate of the frequencies 
over which transport is collision-dominated and the pre- 
ceding calculation is valid: in order for collisions to pro- 
duce relaxation, we require uj ^ t~^, which occurs for 
huj/ksT < Na'^. 

Conductivity with impurities.- We turn now to the con- 
ductivity of the noninteracting, disordered system. We 
restrict to the case of scattering off random point impuri- 
ties, characterized by Vi{r) ~ Vo6{r) and the locations of 
which we shall assume are uncorrelated, {{pi{r)pi{r'))) cx 
6{r — r'). With these assumptions, we are led to 
Hd in © with U{r) = J dr'vi{r ~ r')pi(r'). As be- 
fore, we first compute the conductivity for a single 
node. Assuming that the impurities are sufficiently di- 
lute that the Born approximation is valid, the quasi- 
particle lifetime due to impurity scattering from a sin- 
gle node is given by = — 2ImI]''''*(w, fc) where 



is the re- 



tarded self-energy, G^^\uj, k) — (lo + i6 — Xvpk)"^ is the 
Green's function for a noninteracting Weyl fermion with 
helicity A, and the form factor from the overlap of helicity 
eigenspinors, J'\\'{k,k') = ^{1 + XX' cosOf^i^/) to leading 

7^ = 27r7g(a;), where 



order. We find (see App. IC 1 

g{uj) — 2tx'^v-^ density of states and 7 = ^rixv^ 

characterizes the strength of the impurity potential. 
To evaluate the conductivity we use the Kubo formula. 
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where nj^°*(w, q) is the retarded response function which 
(7) for a system of linear dimension L is defined to be 



n;r.'(^,'7) = - 
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with = {t,r), p^ ^ and = {-ii^-i]:.!). From 

gauge invariance W^l{ijj,q) = 11'°* (w, |q|) - 2^ 

so that (suppressing q — 0), <j[ijj^T) — — ^Imn''°*(w) 
— ^lmIF^^{ijo) . Some algebra yields (see App. IC 2} 
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where frii^) — [e"/"^ + 1]^^ is the Fermi-Dirac func- 
tion and we have used the retarded helicity-basis Green's 
function dressed with disorder lines, G5f*(i^,A;) = [ui — 
After a tedious calculation, we 

'^J{Lo,f) where f = T/uo, 



Xvpk + i/2T{uj) 
may write cr(a;,r) 
UJ — uj/ujo, so that fri'^) 



ff{uj), ujQ = 27rw|/7 is 
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Figure 1. Frequency-dependent conductivity of a single Weyl 
node with disorder (constants defined in the text.) 
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Figure 2. pdc = cr^^ '^t (defined in the text) for the inset 
parameter values compared to experimental data from^. 



a characteristic scale set by the disorder strength, and 
J{Lo,f) = l/f l:^ii£b|d£±Mi(e + ^,e)_withIaconi- 



plicated rational function shown in App. IC 21 

In our model disorder can scatter between nodes, so 
1/t{uj) acquires a factor of N when iV > 1; in common 
with the interacting case, a also has an overall prefactor 
of N . From these it is easy to show that for N nodes, 
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which is identical to the = 1 result (Fig. [T]) upon 
rescaling ujq ujq/N. 

While in general we integrate (|12p numerically, in cer- 
tain limits an analytic treatment is feasible. For u <^T, 
ff{i) ^ ff{i + Lj) « ^iljf'{e). Expanding X in powers 
of UJ and resumming only terms dominant as e — ^ 0, we 
recover the result of Burkov and Balents (Eq. (15) oft): 
namely, a Drude-like response with a width vanishing as 

NT"^. and a finite dc limit of -rr^. 

3/17 

In the opposite limit, T — )> at finite uj we may replace 
the Fermi functions by step functions, which yields 



(uj) « A^ 

^ ' 12hvF 
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The leading term is universal and independent of disor- 
der, and is simply (Jq^\ Both regimes are captured in 
Fig. [1] which shows a{uj,T) for w < beyond which 
the Born approximation is insufficient. 
Experiments.- Ini^ the dc resistivity of polycrystalline 
¥2^267 was found to vary with temperature as pdcT ~ 
130n • cm • K over lOK <T< 170K, which is reminis- 
cent of our result with interactions ([T]) . Accordingly, we 
compare this data with a model of a clean WSM with 
A^ = 24^, as shown in Fig. [2l We find rather good agree- 
ment with experimental data for physically reasonable 
parameter choices, shown inset. Very recently, transport 
in single crystals of another pyrochlore iridate, Eu2lr207, 
has been studied^^ under pressure for 2K < T < 300K; at 
low pressures'^ 2.06 — 6.06 GPa, Pdc{T) resembles Fig. [21 
consistent with WSM behavior. 



Conclusions. - The conductivity of WSMs thus exhibits a 
rich variety of behavior on varying frequency and tem- 
perature, in both the interacting clean and noninteract- 
ing disordered limits, as shown in Figs. [T] and [H In 
particular, its nontrivial dependence on A^ is sensitive to 
the strength of the interactions; with just disorder, we 
find a striking difference between the u <^T and uj ^ T 
regimes, with the T — > ac response dominated by a 
universal, disorder-independent contribution. While the 
limited existing dc conductivity data on the candidate iri- 
dates broadly agrees with our theory in the clean limit, 
we caution that more dc and ac conductivity measure- 
ments on single crystals with controlled disorder are re- 
quired to make a rigorous comparison. Theoretically, the 
interplay of disorder and interactions, and corrections to 
the isotropic node approximation still need to be consid- 
ered. In particular, it would be striking if the distinct 
behavior of the disordered system across the different 
frequency regimes survives the inclusion of interactions; 
these questions we leave open for the future. 
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Appendix A: The Collision Operator and the 
'Leading Log' Approximation 

In order to define the collision operator C, introduced 
in the main text, we need to determine the tree-level 
scattering amplitudes due to the interaction and how 
they couple the quasiparticle distributions at different 
momenta. An additional complication stems from the 
fact that of the A^ Weyl nodes, N /2 have chirality —1. It 
is easily verified that the single-particle dispersion of the 
a*'' Weyl node is given by E\^a{k) = Xa^vpk, so that for 
Xa = the eigenspinors of Ha with eigenvalues ivpk 
are interchanged relative to Xa = +1- For convenience, 
we define Ca = XaA, so that the energy is Ea{k) = eaVpk 
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with = ±1. As the Coulomb interaction is instanta- 
neous, the cohision operator is diagonal in frequency. Hi 
is then given by 

Hi = \l Tx»J4,,(fcife2g)x 

^ J kik2q 

llt,{ki + 9)7^(^2 - q)l.,a{k2)leAki) (Al) 

where 

[Ul{k,-q)U^Sk2)]^^^^ (A2) 

is the scattering amplitude for two-particle scattering 
from {(&, ei, fci), (a, £2, ^2)} to {(a, 63, ^2 - q), {b, £4, fci + 

I 



q)} and in addition depends on the chiralities of the ini- 
tial and final nodes. This dependence comes from the 
fact that the scattering amplitude depends explicitly on 
the overlap of the initial and final spinors, which in the 
e-representation introduced above depends on their chi- 
ralities. We observe here that in the (7 — > limit (which, 
as we will shortly see, dominates the current relaxation) 
T++ ^ ^ = ^ ^ , corresponding to the fact that q ^ 
scattering is indifferent to the chirality, so the fact that 
half the Weyl nodes have a chirality opposing the one 
under consideration is unimportant. 



Using Fermi's golden rule, it is now straightforward to 
show that the collision operator is 



g(k, uj)k 
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j \^5{k-k' -\k + q\ + \k' - q\) f_{k)fl{k')fl{\k + q\)f_{\k' - q\)R, (fc, fc', q) X 
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6{k + k'-\k + q\-\k' - q\) fl{k)f_{k')fl{\k + q\)fl{\k' - q|)i?2 (fc, fc', q) 
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(A3) 



where 



{k,k',q) = \T+^_4k,k\q)-T++^Jk,k',k' -k-q)\ 
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|T++_+(fc, k', q)f + |r|++_(fc, k',k'-k-q)\ 



N 



|r+r_+(fe, k', q)\' + |T+r+_ (fc, fc', k'-k-q) 



(A4) 



and 



i?2 (fe, fe', g) - i |r++++(fe, fc', q) - r++++(fc, fe', k' -k-qf+(^^-?j |T++++(fc, fe', q)|'+ (^^^ |r|+++(fc, fc', g) 

(A5) 
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describe scattering between oppositely and identically 
charged particles, respectively. The diagrams corre- 
sponding to scattering processes are the same as in 
Ref. M- 



Turning now to the C-dependent portion of the expres- 
sion for the functional Q, we may use the symmetries of 
C under combinatorial exchanges of momenta to write (in 
rescaled variables k/T,Ld/T) the first term in Q as 
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/ = / {S{k-k' ~\k + q\ + \k' -q\) fUk)fl{k')fli\k + q\)f_i\k' -q\)R,{k,k\q) 

Jk " Jk.k'q 



^■g{k,u:) - ^g{k',uj) - "|" ^| g(|fc + g| , a;) 



k' 



T9iW - q\ 



S{k + k' -\k + q\- \k' - q\) fjk)f_{k')fl{\k + q\).fli\k' - q\)B.-2 (fc, k\ q) 
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We now specialize to a perturbation of the form 
g{k,uj) = fc^(cj), which is a zero mode of the i?2 pro- 
cesses. This leaves only the terms in (IA6[) proportional to 
We may represent the energy-conserving (5- function 

as 



S{k-\k + q\+n)S (k' - \k' -q\+n) dfl 



Next, we choose coordinates in which q is along z so that 
> (^k - v^fc2 + g2 + 2kqcos9 + iij 

-!^s{cose-^^^±^)eik + n) 



6{k-\k + q\ + n) 



and similarly for k'. The angular integrals over 9 and 
9' can now be done trivially, at the cost of leaving a 
complicated dependence on fl in fc, k' . For q < k, k' , 
the 9 and 9' integrals give one if —q < fl < q and 
otherwise; we will show that it suffices to consider small 
q, so the knowledge of this regime is sufficient. With 
these observations, we have 



gCg 
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dn'^^f_{k)fl{k')fl{\k + q\)f_{\k' - q\)R, (fc, k', q) 

k.k' ,q J -q "-9 '^9 



(A7) 



The cross-helicity scattering rate Ri contains terms of 
the form V{qY, which diverge as 1/q^ at small q for di- 
rect scattering; the exchange contributions remain finite 
and do not contribute at leading order. The integral over 
the dummy 17 variable gives one power of q in thenumer- 
ator and substitution of our ansatz shows that the term 
in braces gives two more powers of g, so that at small 5, 
the integrand is proportional to 1/q. which gives a log- 
arithmic divergence. It is this logarithmically dominant 
contribution that is referred to in the text as Cq, whence 
we implicitly define Ci = C — Cq. 

Having identified the dominant contribution to the col- 
lision integral, we proceed as follows: first, we substitute 
our ansatz into the expression above, and perform the in- 
tegral over g; as this integral is logarithmically divergent 
in both the UV and the IR, we will need to regulate it 
appropriately. The UV cutoff is the temperature T, as 
this is the largest physical energy scale in the problem; in 
the IR, the divergence is avoided due to thermal screen- 
ing of the Coulomb interaction which occurs at an energy 
scale set by both the coupling constant and the temper- 



ature, which we take to be aT. With these choices, we 
may perform the integral over q. and then set q = in 
the remainder, leaving integrals over $7, k and k' which 
may then be performed straightforwardly; upon includ- 
ing the parts of the functional Q that are independent of 
the collision operator, we arrive at the results quoted in 
the text for the LLA. 



Appendix B: Renormalization of couplings 

As mentioned in the text, in order to make a mean- 
ingful comparision between the calculated conductivity 
and experimental data at temperature T, we need to 
use values of the various parameters renormalized to 
the characteristic energy scale ksT. In order to deter- 
mine these, we work with the Euclidean path integral 
Z[A] = ^[DAoDijj''Dip]e-^^'l''^°'^\ written in terms of 
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Figure 3. One-loop diagrams that govern the flow of the 
marginal couplings vp and e. 



the gauge-fixed action 
S[ij,Ao,A] ^ J drdr 



iPi{r,T)[dr - i\e\Ar{r,T) 



+VF(T ■ (iV + \e\A/c)]ipa{r, r) 



47r 
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which (as is easily verified) is equivalent to the Hamilto- 
nian H described earlier. Here, A is a static external vec- 
tor potential, while A^- is a dynamical field that mediates 
the /{er) instantaneous Coulomb interaction between 
the electrons. We implement a momentum-space cutoff 
< A, but leave the frequencies unrestricted so that the 
interaction always remains instantaneous; we then study 
the flow of S as we lower A. The action has three dimen- 
sionless parameters - the Fermi velocity wf, the electric 
charge e and the dielectric constant e - all of which are 
unrenormalized, and hence marginal, at tree level. Of 
these, e is further protected against renormalization to all 
orders in perturbation theory by gauge invariance. How- 
ever, vf and e have no such protection and at one loop 
they acquire corrections of 0{a^) when A is lowered 
to A — (5A. The diagrams governing the flow of vp and 
e at this order in a are the one-loop self-energy and the 
fermion (RPA) bubble (see Fig. [3]) diagrams respectively; 
performing the sum over flavors where appropriate, we 
find that the corrections are, respectively, independent 
of and proportional to N . Evaluation of the correspond- 
ing integrals leads to the flow equations 



dvF 



2 de 
avF 



N 



da 



di Stt 



as 



N + 2 
Stt 



where £ — log A is the usual RG scale parameter. In- 
tegrating (|B2p between an upper cutoff (which we leave 
unspecified for the time being, but will take to be much 
higher than any scale relevant to transport measure- 
ments, but lower than the energy at which band cur- 
vature corrections become significant) and the tempera- 
ture scale, T gives the renormalized parameters used in 
the text. These renormalized values oi vp and a are to 
be used in the formula for CTdc [Eq-(l) of the main text] 
to compute the conductivity at temperature T, while e 
retains its bare value due to gauge invariance. Represen- 
tative results are shown in Fig. [5] 
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Figure 4. Resistivity as a function of temperature for A'^-node 
Weyl semimetals with A'^ = 1 and A'' = 24. 
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Figure 5. Renormalized fine-structure constant as a function 
of temperature for A-node Weyl semimetals with A' = 1 and 
N = 24. 



We remark on an important distinction between the 
case of the Weyl semimetal in d = 3 and its two- 
dimensional analog, graphene. In the case of graphene, 
while the electrons are restricted to two dimensions, their 
interaction continues to be of the 1 jr form since the elec- 
tromagnetic field lines are not thus constrained; in mo- 
mentum space, this leads to a non-analytic form for the 
instantaneous interaction, V^(q) ~ l/l?!- Since the RG 
procedure outlined above can only produce analytic cor- 
rections at any order, e thus remains unrenormalized. In 
d = 3, we have an analytic ^(q) ^ l/'?'^, and so e can 
and does receive corrections at one loop and beyond. 



Appendix C: Conductivity with Disorder 

1. Scattering Time 

The scattering time is given in terms of the imag 
inary part of the retarded self-energy by^:^ = 



8 



-2ImS5f*(w, fe), where 



we arrive at the result, 



(CI) 

with all quantities defined as in the text. Imposing a 
cutoff |fe| < fl/vp, performing the angular integration, 
and rewriting the negative-helicity Green's function by 
changing variables we have 



de [hie]- hie + 
27r Lo 
(Pk 



, r Ml. 

^ / 7233li^Gr(e + c.,fc)ImGl<=,*(6,fc) 

A, A"' ^ ' 



(C6) 



ImI]5^'=*(cj,A;) = 



Im 



? ? 
riiVQUj 



independent of A, whence 



dx ■ 



iS 



(C2) 



2Tr^g{uj), where g{uj) 

,,2 



2^2yi is the density of states and 7 = |?^i^'o• Note that 

Rel]''''*(aj, fc) is formally divergent with the cutoff regu- 
larization used here; however this can be absorbed into a 
renormalization of vp and we shall not discuss it further. 



In order to perform the integrations in (|C6I) . it is con- 
snient to wri 
tion, we have 



venient to write „ ] , = ^ 1 oj'^ = —. With this defini- 

2r{uj) 2iTvi. uja 



ImG^°'(w,A;) = Im 



1 



uj — Xvpk + i- 



1 



(C7) 



2. Evaluation of the Kubo Response 



In order to evaluate the conductivity we require the re- 
tarded response function nj]'J(cj) which is obtained most 
simply by working in the Matsubara formalism and an- 
alytically continuing the results to real frequency. The 
leading contribution to the Matsubara response func- 
tion corresponds to the particle-hole bubble with 
impurity-dressed propagators, and in terms of Matsub- 
ara Green's functions in the pseudospin basis takes the 
form 

npi,(icj„) = — ^ ^ J [Giii^^n + Cm), 



/5 Jk 

X Q{iem,k)at,] 



(03) 



where Q {iuim k) = (iuJn — <^'k+i/2T{uin))^^ and the trace 
is over the pseudospin indices. Performing this trace and 
summing over the spatial index fi, we obtain 



/ ^?a(«(w„ + er„),fc)^A'(iem,fc) 



A, A' 



(C4) 



where QxiitOni k) — {iujn ~ Xvpk + i/2T(w„))~^ is the 
Matsubara Green's function, now written in the helicity 
basis. Performing the sum over Matsubara frequencies 
by standard methods, analytically continuing the result 
to real frequencies via LOn ^ uj + iS and using 



cr(tJ,T) 



1 



UJ 



Imn''<=*(a;) 



1 

3aj 



imn-;(^). 



(C5) 



By a straightforward change of variable, followed by an 
angular integral we can write the sum over helicities and 
the dimensionless frequency parameter cD^ — uJi/ujQ, for 
convenience replacing e + uj and e by uJi 



E 



lmGl'^iui,k)lmG\^^iLj2,k) 



27r7 



I{UJI,UJ2) 

(C8) 

which serves to define the function X, which contains 
a single integral over the dimensionless parameter x — 
vp\k\/ujQ. Rewriting the Kubo formula and reinstating 
factors of H by dimensional analysis we fin we can express 
the conductivity in terms of a single function of rescaled 
temperature and frequency. 



cr(w,T) 



/17 



'-JiuJ,f) 



(C9) 



where T = T/ujq, uj — uj/ujq, so that fx^uj) — fj,{uj), 
ujQ — 2Trvp /j is a characteristic scale set by the strength 
of the disorder, and 



J(uJ-T) 



de [ft{e)-ff{i + uj) 



+ uJ,e) 



2tt uj 

(CIO) 

Numerical integration of this expression yields the con- 
ductivity for all T, cj, as shown in Figure 1 of the main 
text, and direct evaluation of the integrals yields an ex- 
plicit form for I, 
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X(e + CJ, e) = 



(Cll) 
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We have used the fact that at particle-hole symmetry, the 
deviation from equilibrium is proportional to the sign of A, 
and/0(fc) = l-/^(fc). 

As mentioned in the introduction, a and vf take on their 
renormalized values at temperature T. 



